Report for Model 2

— for both jointly and separately

2024-08-19

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.00       1.00
## dh0           1.76       3.14
## dh1           2.08       4.06
## dh2           1.22       1.76
## dh3           1.01       1.05
## dl0           1.05       1.11
## dl1           1.44       2.32
## dl2           1.10       1.38
## dl3           1.01       1.02
## 
## Multivariate psrf
## 
## 1.9

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1218.56835 22055.6116
DIC3 1159.04687 22132.8959
PWAIC 43.50255 239.8298
WAIC 1186.08453 22153.7768

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.00       1.01
## dh0           1.01       1.03
## dh1           2.20       6.96
## dh2           1.21       1.75
## dh3           1.59       2.98
## dl0           1.02       1.07
## dl1           1.08       1.33
## dl2           1.01       1.03
## dl3           1.04       1.16
## 
## Multivariate psrf
## 
## 1.95

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1380.31317 22187.433
DIC3 1254.95288 22235.386
PWAIC 93.24739 295.985
WAIC 1325.63251 22267.482